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SUMMARY 



Design of a control law for simultaneously suppressing the symmetric and 

0f a Sdng “ model 

IS described. The flutter suppression control law was designed using linear auadrarir 

Gaussian theory, and involved control law order reduction, a gain root-locCs study and use 

of previous experimental results. A 23% increase in the op^n-loop flutterXS Pressure 

syr^etri^ ^‘tunnel test Rapid roll maneuvers at 11% abovedie 

configuration b ° anduy Were aso P erformed when the model was in a free-to-roll 


INTRODUCTION 


A summary of the Active Flexible Wing (AFW) Program is presented in Ref 1 w;rain 
^,?rd a A n ^ ange °f thC Lang,e y Research Center' T ra S 

fixed-in-roll configured™ m0 H d h&d and antisymmetric flutter modes, in a 

freSo roH confiaS^T ^ & ^ etnc flutter ^ only, when the model was in a 
. to ™ configuration. The active flutter suppression system (FSS) test goals were to 

demonstrate: a) simultaneous symmetric and antisymmetric flutter suppressfon for the 

co n t^ta C Z fi Sr™f nd ? flUtKr -ppressionVffiSmn 

ZSI FSS^hSL?f W f '° if 5 ' a r ? ling maneuver '<*><1 alleviation system 
fliitter^HtSfr ^ h 0pen ' loo P flutter boundary. Since the free-to-roll symmetric 

characteristics ^a sin g^FS S rrCufCC ^ symmetric flutter had very similar 

cnaracteristics a single FSS control law was designed and demonstrated for both the flutter 

test configurations, a) and b) as stated above. This paper addresses die' madiematicaf 
modeling, control law design and wind-tunnel test results mathematical 


NOMENCLATURE 


A, B 
C, D 


Bo 

Co 

E 

F, G 
G w 

g 

H 


control law state-space matrices 
control law output matrices 
Kalman state estimator gain matrix 
optimal regulator gain matrix 
expectation operator 
plant state-space matrices 
gust input matrix 
gravitational acceleration constant 
sensor output matrix 
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IN? i 



I identity matrix 

M Mach number 

p estimator Riccati equation solution 

q dynamic pressure, psf 

q f flutter dynamic pressure, psf 

Ql plant output weighting matrix 

Q2 control input weighting matrix 

R v measurement noise intensity matrix 

r w gust input noise intensity 

S regulator Riccati equation solution 

s Laplace variable 

T sample period, seconds 

t time, seconds 

u control input vector 

v measurement noise vector 

w gust input noise 

(0 frequency, radians/second 

x plant state vector 

Xc control law state vector 

y measurement vector 

z accelerometer output, g's 

6 control surface angular position, degrees 

Subscripts: . . . 

LEI wing leading edge inboard 
LEO wing leading edge outboard 
TEI wing trailing edge inboard 
TEO wing trailing edge outboard 
tip wing tip 
Abbreviations! 

AFW active flexible wing 
CL closed loop 

CPE controller performance evaluation 
FSS flutter suppression system 
LQG linear quadratic Gaussian 
OL open loop 

psf pounds per square foot 
RMLA rolling maneuver load alleviation 
rms root mean square 
SISO single-input single-output 


AFW EQUATIONS OF MOTION 

and 
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sssszr tam a fmi,e - 

State-space Equations 

A set of state-space mathematical models were developed 3 for control law d™™ tw 
four S a ™^c7ae 0 |^s e f™ 0 ^rn a, ^S osc “ at< ?> a ^y"uniics approximation usrt 
included cotrecttons for control surface effectiveness based TOS ^ind. 


and 


dx/dt — Fx + Gu + G w w 
y = Hx+ v 


( 1 ) 

( 2 ) 


SSr«! s; r !!tt s« 

^tR's.TSsassr*--’* 

Open-loop Dynamic Pressure Root-locus 

200 250^ models at six dynamic pressures, q = 100, 150 

stewed' T^n , 3 i? n P H f J thC flexible ' m °de root-loci with dynamic pressure were 

1 he open-loop, dynamic pressure root-locus of the first four flexible symmetric 

res^S'S^ “VS? f “ ed -' n : ro " configuration, are shown fo fi^Tnd 3 

Sg ,he ""““-jnstabiiity. The unstable mode wafprimarily 

FS ° r ll; fc> L ooth the symmetric and the antisymmetric motions. The sixth and 

2) At Mach”? 5 the maUniS* 16 #t | Ue “ ci ' S als0 Knded to coalesce (not shown in figure 
esti™!Si m £ 248 off a n 1 H7 n Tl°° P S ?“ flu ,' “dynamic pressure was 

pressure wasestirafed to bills prfat f 09 Ite XcS flutter d y namic 

locus is also shown in figures Ian'S aid ^—-Pressure root- 


CONTROL LAW DESIGN 


laws™ e hSrSS!«nS" ° bjec,ive 10 devel °P '°»-°foer robust digital control 
of the rnndel in t£r y s “PP ress the symmetric and antisymmetric flutter modes 

of the model in the fixed-in-roll configuration with allowable control surface acdWtv^Te 

maximum pemussible control surface rms deflection and rates were 1 .0 degree (IITl ?Hz 

f U ter frequenc y) and 75 degrees/second, respectively. From the 1989 test 2 the 

conSrr 010 flU f ter fr6qUenCy was known to be 1 .8 Hz below ^fanalytical value. The 
w was also required to be sufficiently robust to compensate for this difference. 
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2 a t, 3) ' 

la^^is slrown 

nassed through 25 Hz first-order antialiasing ? f ^fyi Hz Thedicital 

i£S^ 

and left actuators after 0.005 seconds computational delay. 


Design Plant Model 

The 68th order antisymmetric state- space^uation ^ anLysi^and^he 1989 

configuration was used as the design P ^ rn e » si rilica i anc [ was encountered at a 

^®SS5S3^.-as5« 

ESSS ESS 

as candidates for measurement inputs and control outputs, respecti y. 


(3) 

(4) 


Full-order LQG Design 

equ ^<^ 

(4), minimizes a weighted quadratic cost function defined by E[y Qiy ^2 1. 
and Q2 are the plant output and control input we.ght.ng tnatnees • . 

dxjdt = A 0 x c + B 0 y . 
u = Cox c , 

A 0 = [F-BoH + GCo] 

B 0 = PHTRv' 1 
C 0 = -Q2- 1 G t S. 

steady state dual matrix Riccati equations, given by 


where 


FP + PFT + G w RwGw T -PHTR v -'HP = 0 

SF + FTS+HTQ.H - SGQj-'GTS = 0, 
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where R w and R v denote the intensity matrices of the gust input and me^nrempnf 

M “rfer cJnLnaw ' and ° f and ^sity lattices for®, he 

■ d the subsa J uent order reduction process were determined after 

S ^"T 8 '7 0rder COn ^ lkr ™ foand“ot„al 

8 pBUH moael - lhe contro1 law order reduction process is described next. 


uraer Reduction 

Q j //^\ » l (4), in block-diagonalized form, are shown in erniatinnQ 

and (6), where .he vector xc represents the retained states and tta?«Tr xc 2 regents the 
remaining states assoc, ated with die damped higher frequency dynamics 


d_ 

dt 


K,i 

. = ' 

l x ft 



A , o 

O 01 A 


o2 


' X cl 
X c2 + 


,ol 

o2 


U = 


[ C ol C o2 


ft 


(5) 

( 6 ) 


dynamics “n^ 

denvanve dx^d, to zero and solving for x c2 , provided the matrix Ao 2 "f nonI!nllar4. 
uce state space model of the control law is given by equations (7) and (8). 

dxc/dt = Axc + By 
u = C x c + d y 


where 


(7) 

( 8 ) 


x c = x c 1 , B = B 0 1 , C = C 0 1 
and D = -C o2 Ao2- 1 Bo2. 


i^aaaaaaaasar* 

ssSSSSSSSSSSH 

’ wo ' in P ut two-output control law, is given by equations (9) and (10). 


dx c 

dt 



-5.2 64.6 

-64.6 -5.2 

-0.4 2.1 
3.6 -9.4 


x c + 


x c + 


1.25 1.95 
-0.45 -0.73 


-0.06 

0.13 


-0.09 

0.21 



(9) 


( 10 ) 
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The cotresponding Bode diagrams of the 

H h TrlT<ffi& a^ iraTfo ”te V primary stabilizing gain of Uris control law 
Jas from the sensor > to the control surface 

the symmetric and antisymmemc plant models at 3 W ps^tesKp ^pons ^ 
eliminated. However, with the addition ofT ^.TO5 wcond compn^ m i ; nally stabte . 

(SlSO)conoVlwb, slsoiCa ia" taiha. 

fflSSlSSKS compmate for the Lnputatmnal delay effects, and possible 
un^in^nTe actual fiuner frequencies, as mennoned earher. 


SISO Control Law 

. cimnlified SISO control law (plot labeled by SjEofaipa in figure 5) was improved 

half plane, with a negative feedback gain of 1.3 level 

actuator poles near 50 Hz bwome unstable with subsequent gain attenuation at 

of at least 1.3 deg/g in the 8 to 12 Hz eque y S y St em, and accommodate the possible 
higher frequencies, was necessary to sta . flutter freauencies In addition, 

difference between the filter and one cycle computational 

£ £» these two effects was about 40 

degrees at the frequency 10 Hz. 

The gain and phase compensates were 
and D in the SISO control law, and studying th gmn and phase diagreuns 
loop stability responses. An increase ir .Cl i and ^asem resu^ ^ ^ ^ 

increase at low frequencies. An increase oarameters were varied, until a gain- 

frequencies, which wasalsobeneficial. The P . nge 8 to 12 Hz, and 

level near 1.3 deg/g (2.3 dB) law complex pole was also 

sufficient phase lead was obtained. P The high frequency gain was 

ke^ro^ ^ 

(12), assuming negative feedback. 


dx 


C __ 


dt 


-6.0 64.6 
—64.6 — 6.0 


x c + 


1.95 

-0.73 


tip 


(ID 


8 = [ 14.4 -3.1] x c + 0.63 z, 

°teo L 


up 


( 12 ) 


140 



nziln* corresponding gain and phase plots are shown in figure 5 and are labeled &rEo^tip 
(MbO). The complex poles and zeros of this control law were -6±j64.6 and -30+j56, 
respectively. A second-order notch filter, given by the transfer function (s 2 +42s+44100)/ 

u + ^. + ^. 100 was added t0 increase the symmetric model gain margin to 6 dB, near 33 
Hz. Tnis filter attenuated a 33 Hz lightly damped oscillation due to the interaction of the 
sixth and seventh symmetric flexible modes. A first-order washout filter, given by the 

signal funCtlon s ^ s+6 ^’ was aIso added t0 remove any steady state input bias to the sensor 


Discretization. The resulting 5th order SISO control law in Laplace domain was 
discretiz^ using the Tustin transformation z = (l+sT/2) /(l-sT/2), where T is the sampling 
the 200 Hz sampling rate used by the digital controller, T = 0.005 seconds 
With the Tustin transformation at this sampling rate, the Bode diagrams in the Laplace 
domain and the discrete domain were almost identical below 15 Hz. Hence no frequency 
warping corrections were applied. 


Dynamic-pressure root-locus: The open- and closed-loop dynamic pressure root-locus 
plots are compared in figures 2 and 3. These comparisons indicated that both the symmetric 
and antisymmetric models were stable, up to dynamic pressure q = 350 psf. The closed- 
u° P T^ eC l UenCy decou P lin g wa s due to lowering of the frequency of mode 2 to about 6 8 
Hz. The frequency of mode 3 was increased to 1 1.6 Hz, but the damping ratio was oniy of 
the order 0.010 at 300 psf. 


Sensitivity studies. The closed-loop system sensitivity was studied by perturbing the 
SC i^ d j l d modal frequencies in the state-space block-diagonalized plant model by 
±10% and the nominal gains by ±4 dB at q = 250 psf and examining the closed- loop 
system step responses, for all possible combinations. These studies indicated that the 
design could accommodate simultaneous gain and frequency changes for all cases except 
when the second and third mode frequencies were perturbed to approach each other. 
Sensitivity studies were also done using the state-space model with and without the 25 Hz 
antialiasing filters, with and without one cycle delay, with additional delays, and with ± 6 
do gain perturbations at 250 psf. These studies indicated that the symmetric configuration 
could tolerate one additional delay (or phase lag of 1.8 degrees/Hz) at half the nominal 
giun, but the antisymmetric configuration would become unstable with an 1 1 Hz 
oscillation. The phase and gain margin comparisons with the experimental results, 
described in the next section, indicated that this particular situation may have been 
dia^arrf^ ex P er ’ ment - The gain loss was apparent from the experimental Bode 


SUMMARY OF TEST RESULTS 


Open-loop Flutter. Based on examination of the peak-hold data obtained during the wind 
tunnel test with the tip ballast store coupled, the open-loop (OL) flutter dynamic pressures 
were as follows: The free-to-roll OL symmetric flutter was at a dynamic pressure of 235 
psf, at a frequency of 9.6 Hz. The fixed-in-roll OL antisymmetric flutter was at a dynamic 
pressure of 219 psf, at a frequency of 9.1 Hz. These experimental symmetric and 
antisymmetric OL flutter dynamic pressures were, respectively, 13 and 14 psf below the 
predicted values, and the flutter frequencies were, respectively, 1.6 Hz and 1.8 Hz below 
the predicted values. 
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Open-loop frequency responses. Figures 7 and 8 show the OL frequency responses of 

zap due to 8 teo from analysis and experiment at 250 psf for the symmetnc and 
antisymmetric (fixed-in-roll) cases, respectively. At this dynamic pressu^ the OL plant 
unstable. So, the OL frequency responses were computed from closed-loop (CL) 
experimental data, using the Controller Performance Evaluation (CPE - ) P^cdure. 
Figure 7 indicates good agreement below 9 Hz and qualitative agreement above • 
Above 12 Hz, the magnitudes differ by about 5 dB while the phase angles 
Figure 8 indicates fair agreement, below 7 Hz, and qualitative agreement above 12 . 

AtLve l2 Hz, the magnitudes differ by 6 to 8 dB and the phase angles differ by 10 to 20 
degrees! Note, that fofeach phase diagram, the 180 degree crossing occurs near the 
respective OL flutter frequencies, and the difference between their predicted and 
experimental values is quite apparent. 


Closed-loop Tests 

The active flutter suppression test results are summarized in figures 9 through 13. 
Figures 9 aSd C 10 show the wind-tunnel test dynamic pressures versus the . free stream Mach 
number During the wind-tunnel test, in the fixed-in-roll configuration, with both die 
symmetric antisymmetric FSS control laws operating, the CL system was stable up to 
q Y = 270 psf, at Mach 0.46. This augmented q represents a 23% increase over the OL 

antisymmetric qf. 

During the wind-tunnel test, in the ffee-to-roll configuration, with the qjnmrtric FS S 
control law operating, the CL system was stable up to q = 290 psf, at Mach ?-48- 
augmented q represents a 23% increase over the OL symmetric qf as shown m figure . 
This FSS control law also suppressed the flutter when a Rolling Maneuver Lo 
Alleviation (RMLA 8 ) system was tested with rapid roll maneuvers at q-2wpsf, 0 
above the OL symmetric flutter boundary. This RMLA control law used LEO and 
control surfaces, so the interaction with the FSS control law was minimal. 

The rms deflection and deflection rate of the right and left side TEO control surface 
were computed from the data sampled at 200 Hz at each nxed-iu-rol FSS 
the value of the right and left differed, the maximum is plotted in figure 1 1. The maximum 
rms deflection and rates were less than 0.4 degrees and 25 degrees/second, r^pectively. 
These maximum rms deflection and rate demands of the actuators were well below the 
mSmum Sowable values of 1 deg and 75 deg/sec as stated earlier in the paper. 

The Nyquist-diagram-based gain- and phase-margins were estimated using the CPE 
techniquefduring thfexperiment. These estimates were compared with corresponding 
analytical quantities in figures 12 and 13, for the symmetric free-to-roll and t e 
antisymmetric fixed-in-roll configurations, respectively. For the symmetnc fr^-to-roll 
configuration (figure 12), the analytical and experimental gam margins were above ±6 dB 
ud to 270 psf The analytical positive phase margins (at or below 7 Hz) were about 20 
degrees, but the negative phase margins (at or above 12 Hz) were ^ ell ^ V | 7 Q 5 J gr 
The analytical phase margins were close to experimental results up to about 270 psf. 

For the antisymmetric, fixed-in-roll configuration (figure 13 ), the analytical neg^ auve 
gain margins were only -3 dB.The analytical positive phase m^msfat orbelow 7 Hz) 
were about 20 degrees, but the negative phase margins (at or above Hz) were 5 
degrees. The analytical phase margins were close to the expenmental data at 25 p , ^ 
hecause the design model was fairly accurate at frequencies below 7 Hz (see figu e ). 
negative gain aid phase margins a! the high frequency end were primarily responsible for 
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preserving the system stability. The source of additional phase lag with increasing dynamic 
pressure was possibly due to highly loaded actuators. The gain loss was apparent from the 
experimental Bode diagram shown in figure 8 in the 8 to 12 Hz frequency range. 


CONCLUDING REMARKS 


A single-input single-output control law was designed for flutter suppression using 
linear quadratic Gaussian theory and involved control law order reduction, a gain root- 
locus study and use of previous experimental results. The control law was digitally 
implemented and tested. Simultaneous suppression of symmetric and antisymmetric flutter 
modes m close proximity was demonstrated to 23% above the open-loop antisymmetric 
utter boundary when the model was in a fixed-in-roll configuration. Symmetric flutter 
suppression system operating simultaneously with a rolling maneuver load alleviation 
system was tested to 23% above the open-loop symmetric flutter boundary, when the 
model was in a free-to-roll configuration. With this combined system, rapid roll maneuvers 
were also performed at 1 1% above the symmetric flutter boundary. 
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2 ( 


rigurc 1. Accelerometer and control surface locations on AFW 
wing plan form* 


imaginary part 



Figure 2. Symmetric open- and closed-loop dynamic pressure 
root-locus at M=0.5 (arrows indicate increasing 
dynamic pressure). 


imaginary part 



real part 


Figure 3. Antisymmetric fixed-in -roll open- and closed-loop 
dynamic pressure root-locus at M=0.5 (arrows 
indicate increasing dynamic pressure). 



- left 


Figure 4. Digital FSS control law implementation block 
diagram. 



Frequency, Hz 


Figure 5. Bode diagram of reduced, second-order control laws . 
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imaginary pari 



Figure 6. Gain root-locus plot for negative feedback from ztip » 
&tfo at 350 psf, antisymmetric fixed-in-roll 
configuration (x = poles, o = zeros, * indicates gain 
increment by 0.1). 



Phase, 

degrees 



Frequency, Hz 


Figure 7. Comparison of ztip / &TEO Bode diagrams at 250 psf, 
symmetric configuration. 
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Figure 8. Comparison of ztip / &teo Bode diagram at 250 psf 
antisymmetric fixed-in-roll configuration. 



Mach number 


Figure 9. Summary of results for fixed-in-roll FSS wind-tunnel 
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• ■ sensor and control surface for FSS 
Q C3 sensor and control surface for RMLA 



Figure 10. Summary of results for free-to-roll RMLA/FSS 
wind-tunnel test. 
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Figure 1 1 . Maximum 6 teo control surface deflection and rates 
demands for simultaneous symmetric and 
antisymmetric flutter suppression tests. 








